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Abstract. Convincing evidence for long-term variations in the emission properties of the anomalous X-ray pulsar 
IRXS J170849. 0^00910 has been gathered in the last few years. In particular, and following the pulsar glitches of 1999 
and 2001, XMM-Newton witnessed in 2003 a decline of the X-ray flux accompanied by a definite spectral softening. This sug- 
gested the existence of a correlation between the luminosity and the spectral hardness in this source, similar to that seen in 
the soft 7-repeater SGR 1806-20. Here we report on new Chandra and Swift observations of IRXS J170849.0-4009I0 per- 
formed in 2004 and 2005, respectively. These observations confirm and strengthen the proposed correlation. The trend appears 
to have now reversed; the flux increased and the spectrum is now harder. The consequences of these observations for the twisted 
magnetosphere scenario for anomalous X-ray pulsars are briefly discussed. 
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1. Introduction 

The Anomalous X-ray Pulsars (AXPs) are a small group of 
sources which stand apart from other known classes of X-ray 
pulsars. In particular, they all rotate with spin periods clustered 
in a very narrow range (P ^ 5 — 12 s), they have large pe- 
riod derivatives (P ~ lO^^'^ — 10^^"ss~^), and, except in 
one case (Camilo et al. 2006), deep searches for radio pulsa- 
tions gave so far always negative results (Burgay et al. 2006). 
Another important characteristic which motivated the "anoma- 
lous" label (Mereghetti & Stella 1995; van Paradijs, Taam & 
van den Heuvel 1995), is their relatively high X-ray lumi- 
nosity (~ 10'^'' — 10"^^ erg s~^), which cannot be accounted 
for by rotational energy losses alone, and no convincing evi- 
dence for a companion star was discovered so far for any of 
them. These considerations quite naturally led to the idea that a 
non-standard energy production mechanism is involved in their 
emission. 

Many different models have been suggested all along for 
AXPs, such as they are accreting from a fossil disk, formed 
by the debris of the supernova event, or from a very low-mass 
companion (e.g. Mereghetti & Stella 1995; Mereghetti et al. 
1998; Chatterjee et al. 2000; Perna et al. 2000; Alpar 2001). 
On the other hand, many observational properties support the 
idea of these sources being magnetars, i.e. isolated neutron 
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stars powered by the decay of their huge magnetic fields 
(B - IQi^ - IQis G; Duncan & Thompson 1992; Thompson & 
Duncan 1993, 1995, 1996). In fact, if the large observed spin- 
down is interpreted in terms of magneto-dipolar losses, all the 
AXPs seem to have magnetic fields in excess of the quantum 
critical field (B > 4.4 x lO^^G). If this is the case, AXPs 
should be related to the Soft 7— ray Repeaters (SGRs), another 
class of X-ray sources thought to involve strongly magnetic 
neutron stars (see Woods & Thompson 2004 for a recent re- 
view on SGRs/AXPs). In recent years, intense monitoring pro- 
grams revealed several common features between AXPs and 
SGRs, (i.e. short bursts, weak IR counterparts, high energy 
tails; Gavriil, Kaspi & Woods 2002; Kaspi et al. 2003; Israel 
et al. 2003; Kuiper, Hermsen & Mendez 2004), strengthening 
the idea of an underlying relation between these two classes of 
sources. 

AXPs' spectra in the X-ray range are well described by an 
empirical model, made by an absorbed black body (kT ^ 0.3 — 
0.6 keV) plus a relatively steep power law with (photon index 
p 2 — 4), and a hard X-ray power-law tail with P ~ 1. Until 
a few years ago AXPs were commonly believed to be steady 
X-ray emitters (even if hints for variability were akeady found, 
see Iwasawa et al. 1992; Baykal & Swank 1996; Oosterbroek 
et al. 1998) but recently flux changes and spectral variability 
were detected, both long-term and with spin phase (Kaspi et 
al. 2003; Mereghetti et al. 2004; Rea et al. 2005). 
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IRXS J170849.0-400910is a prototypical AXP, with a pe- 
riod of ^ 11 s (Sugizaki et al. 1997; Israel et al. 1999), a spin- 
down rate of'^ 2xl0~^^ss~^, and a soft spectrum (Israel et al. 
2001). A phase-coherent timing solution, inferred thanks to the 
long Rossi-XTE monitoring of this source, led to the discovery 
of two glitches in the last few years, with very different post- 
gUtch behavior (Kaspi, Lackey & Chakrabarty 2000; Dall'Osso 
et al. 2003; Kaspi & Gavriil 2003). In a very recent paper, 
Rea et al. (2005) showed that both the flux and spectral hard- 
ness reached a maximum level close to the two glitches that 
the source experienced in 1999 and 2001, and then decreased 
again in close correlation. Moreover, a long observation taken 
by BeppoSAX during the recovery from the second, more dra- 
matic, glitch revealed evidence for a relatively broad absorption 
line at '~ 8 keV (Rea et al. 2003), not re-detected in more re- 
cent pointings (Rea et al. 2005). The spectral feature has been 
interpreted as a cyclotron resonance feature, yielding an esti- 
mate of the neutron star magnetic field of either 9.2 x 10^^ G 
or 1.6 X 10^^ G, in the case of electron or proton cyclotron 
absorption, respectively (Rea et al. 2003). In this paper we 
present new Swift (calibration) and Chandra observations of 
IRXS J170849.0^00910. In Sect.|2]and[3]we describe the tim- 
ing and spectral analysis of the Swift and Chandra observa- 
tions, respectively, as well as reporting the results. Discussion 
follows in Sect.|4] 

2. Swift observations 

IRXS J170849.0^00910 was observed with the Swift satelHte 
(Gehrels et al. 2004) a few times, being a calibrator for tim- 
ing accuracy and for the wings of the Point Spread Function of 
the X-Ray Telescope (XRT, Burrows et al. 2005; see Table [T] 
for a detailed log of the observations). Here we focus on data 
taken in Window Timing (WT) mode and in Photon Counting 
(PC) mod43 longer than 1 ks. In particular, we used the PC data 
for spectral analysis and the PC + WT data for timing analysis. 
This choice is dictated by the fact that WT observations are af- 
fected by the source being at the edge of the window and dur- 
ing part of the observation the source fell outside the window, 
making difficult a secure evaluation of the instrument spectral 
response. 

Data were analysed with the FTOOL task xrtpipeline 
(version build- 14 under HEADAS 6.0). We applied standard 
screening criteria to the data (CCD temperature T < —45° C, 
eliminated hot pixels and bad aspect times). Hot and flickering 
pixels were removed. High background intervals due to dark 
current enhancements and bright Earth limb were removed too. 
Screened event files where then used to derive light curves and 

' Swift-XKT can observe sources in three observing modes: 
Low Rate Photodiode (LRPD), Window Timing (WT) and Photon 
Counting (PC), having a timing resolution of 0.14 ms, 1.8 ms, 2.5 
s, respectively. Each mode is designed to deal with sources of differ- 
ent intensities in order to minimize the effects of photon pile-up but 
losing spatial information. In LRPD the entire CCD is read as a pho- 
todiode and there is no spatial information. In WT mode a ID image 
is obtained reading data compressing along the central 200 pixels in a 
single row. PC data produce standard 2D images. For more details see 
Hill et al. (2005). 
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a - Observations in bold face are the ones used for spectral and 
timing analysis. Observation in italic face have been used for timing 
analysis only. Observations shorter than 100 s were not listed above. 



spectra. We included data between 0.5 and 10 keV, where PC 
response matrix is calibrated (we used the v.8 response matri- 
ces). 

We extracted data from two WT observations. The extrac- 
tion region is computed automatically by the analysis software 
and is a box 40 pixels along the WT strip, centered on source, 
encompassing ~ 98% of the Point Spread Function in this ob- 
serving mode. We extracted photons from PC data from an 
annular region (3 pixels inner radius, 30 pixel outer radius) 
in order to avoid pile-up contamination. We consider standard 
grades 0-2 in WT and 0-12 in PC modes. Background spectra 
were taken from close-by regions free of sources. 

2.1. Timing Analysis 

Data were barycentered using the FTOOL task barycorr to 
correct the photon aiTival times to the Solar system barycenter. 
A period search led to a clear detection of the neutron star spin 
period. The best period is P = 11.0027 ± 0.0003 s (afl errors 
in the text are given at 90% confidence level). This has been 
derived with phase fitting techniques. This period is consis- 
tent with the extrapolation from known ephemerides at a con- 
stant period derivative (Kaspi & Gavriil 2003; Dall'Osso et al. 
2003). 

We divided the data in four energy bands then folded them 
at the neutron star spin period. The pulse profiles were then fit- 
ted with a sine wave obtaining the pulsed fraction (PF) in differ- 
ent energy ranges. We define here as PF the (semi-)amplitude 
of the best fitting sine to the normalised and background cor- 
rected folded data. We found a PF of 31 ±2%, 39±3%, 29±4% 
and 35 ± 7% in the 0.2-10 keV, 0.2-2 keV, 2-4 keV and 4-10 
keV energy bands, respectively. 
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Table 2. Spectral parameters of the Chandra ACIS-S observation (4.5 count s^^) and of the simultaneous fits of the 5 longer 
Swift observations (0.3 count s^^ on average). Fluxes are in the 0.5-lOkeV band in units of 10^^" erg s~^ cm^^; column density 
were fixed at the XMM-Newton value of 1.36x 10^^ cm"^ (Rea et al. 2005; phabs model in XSPEC); all errors are at 90% 
confidence level. Normalisations are in XSPEC units, i.e. the number of photons kev^^ s^^ cm^'^ at 1 keV. 
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2.2. Spectral analysis 

Spectral modelling was performed fitting together the five PC 
observations with exposure time longer than 1 ks, grouping the 
PC spectra to 60 counts per energy bin. The spectra were fitted 
in the 1-10 keV energy range since the high absorption made 
useless the data below 1 keV (Romano et al. 2005). For the PC 
data described above we generated the appropriate arf files with 
the FTOOL task xrtmkarf and used the latest v.8 response 
matrices. The spectral parameters are reported in Table 2. 

We first fitted all the data with an absorbed (using phabs 
within XSPEC) power-law, leaving all the parameters free to 
vary. This model gave a reduced value of x^cd = 1-08. The 
resulting column density is Nh = (2.00to;^g) x 10^^ cm'^ 
and the power law photon index is F = S.SS^q'^j. We also con- 
sider a fit with the inclusion of a black body component. In this 
case we derive a lower column density Nh — (l-ST^Q ^g) x 
10^^ cm"^ consistent with the value obtained by Rea et al. 
(2005). In addition we got kT = 0.42tJ^ j]^ keV and i? = 5.4± 
1.5 km (calculated at 5 kpc distance) as well as a flatter power 
law index F = 2.73t°oil- 

Also in this case the fit is acceptable 
(Xred — l-Ol)- The inclusion of the black body component is 
significant at 3cr level (based on an F-test). The 0.5-10 keV 
absorbed (unabsorbed) flux is 4.4tQ Qg x 10^^^ erg s^^ cm^^ 
(1.43^0 I5 X 10^^° erg s^^ cm^^), with the power law com- 
ponent accounting for 72'tl% of the total flux. The power law 
component seems to be slightly decreased from the 82 ± 1% for 
XMM-Newton. For comparison with previous spectra we com- 
puted also the same quantities fixing the column density to the 
XMM-Newton \alue(NH,x MM ^ (1.36±0.04) x IO22 cm^^). 
Results are reported in Table|2](in this case the power law con- 
tribution amounts to 71 ± 3%). 

The IRXS 1170849.0^00910 spectrum and flux changed 
significantly from the XMM-Newton observation in August 
2003. Constraining all the spectral parameters to the XMM- 
Newton values within their 90% confidence intervals, the re- 
sulting fit is not acceptable with Xicd — 5-6. In our best fit, 
the blackbody temperature remains consistent with the XMM- 
Newton value. The photon index instead decreased significantly 
and, at the same time, the flux increased. Interestingly, this is 
in good agreement with the correlation found in this source by 
Rea et al. (2005; see also below). 



3. Chandra observation 

IRXS J170849.0-400910 was observed by Chandra on 2004 
July 3 (Obs-ID: 4605), for -30 ks with the Advanced CCD 
Imaging Spectrometer (ACIS). The ACIS CCDs SI, S2, S3, 
S4, 12 and 13 were on during the observation. In order to 
avoid the pile-up, the source was observed in the Continuous 
Clocking (CC) mode (CC33_FAINT; time resolution 2.85 ms). 
The source was positioned in the back-illuminated ACIS-S3 
CCD on the nominal target position. The data were reprocessed 
using CIAO software (version 3.2). A detailed description on 
the analysis procedures, such as extraction regions, corrections 
and filtering applied to the source events and spectra can be 
found in Rea et al. (2005b). 

In order to perform the timing analysis we corrected the 
events arrival times for the bary center of the solar system (with 
the CIAO axbary tool) using the provided ephemeris. We 
used for the timing analysis only the events in the 0.3-8keV 
energy range and the standard Xronos tools (version 5.19). 
One fundamental peak plus one harmonic were present in the 
power-spectrum. A period of 11.00223 ± 0.00005 s has been 
detected referred to MJD 53189. The pulse profile has not 
changed with respect to the previous detection and the 0.3- 
8 keV PF is 35.4 ± 0.6% (see Fig. 2 right panel). 

Being the CC mode not yet spectrally calibrated, the TE 
mode response matrices (rmf) and ancillary files (arf) are gen- 
erally used for the spectral analysis (see Rea et al. 2005b for a 
detailed description of the matrices extraction). 

We fixed the absorption at Nh.xmm = (1-36 ± 0.04) x 
10^^ cm^^ during the spectral fitting, because of the low statis- 
tics below 1 keV and especially because of well known cal- 
ibration issues at 1-2 keV, as previously reported for other 
CC mode observations (Jonker et al 2003; Rea et al. 2005b). 
Actually, to avoid any CC mode calibration problem, all the 
fits were performed removing the data in the 0.9-2 keV range. 

Also for this observation the best fitting model was the ab- 
sorbed power law plus a blackbody. The blackbody tempera- 
ture does not change much with respect to the XMM-Newton 
detection. The black body radius, 3.6 ± 0.4km, is however 
smaller and the decrease appears to be significant at or above 
the 3(7 level. On the other hand the power law contribution in 
this Chandra observation was 80 ± 2%, still consistent with the 
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Fig. 1. Left Panel: Spectra of the five longer Swift PC observations fitted with an absorbed blackbody plus a power-law. See 
Table 2 for further details on the spectral parameters. Right Panel: folded light curves in three energy bands, from top: 0.2-2 keV, 
2^keV and4-10keV. 



XMM-Newton observation of the previous year However, also 
in this case the photon index is moving harder and the flux in- 
creasing toward what Swift saw a year later (see Sec.|2]i. The 
spectral results are reported in Table 2 and Fig. 2 (left panel). 

4. Discussion 

In this paper we present a new Chandra observation of the 
AXP IRXS J170849.0^00910 and the first Swift observations 
of this source performed as a part of the XRT calibration pro- 
gramme. 

We performed spectral and timing analysis of the data. The 
measured periods, even though affected by relatively large er- 
rors (see § 12. Il l, allows us to confirm that the sources is still in 
a phase of steady spin-down, following the last glitch. 

The spectral analysis reveals the source undergoing sig- 
nificant spectral changes. Interestingly, the trend monitored 
following the glitches epochs (Kaspi, Lackey & Chakrabarty 
2000; DairOsso et al. 2003; Kaspi & Gavriil 2003) and until 
the last XMM-Newton observation has now reversed (see Fig. 
O. In particular, the source spectrum became much harder and 
the total unabsorbed flux in the 0.5-10 keV energy band a frac- 
tion ~ 50% higher with respect to that measured by XMM- 
Newton (Rea et al. 2005). Moreover, our analysis indicates that 
the flux increase is mainly due to an increase in the contribution 
of the thermal component, while the power law contribution to 
the total flux shghtly decreased (71 ± 3%, while the XMM- 
Newton measurement was 82 ± 1%). 

In Rea et al. (2005) it has been proposed that the observed 
correlation between the X-ray flux and the spectral hardness 
may be explained within the "twisted magnetosphere" sce- 
nario (Thompson, Lyutikov & Kulkarni 2002; Beloborodov & 
Thompson 2006). The basic idea is that when a static twist is 
implanted, currents flow into the magnetosphere. As the twist 
angle A0ns grows, charge carriers (electrons and ions) provide 
an increasing optical depth to resonant cyclotron scattering and 
hence a flatter power-law. At the same time, the larger returning 



currents heat the star surface producing more thermal photons. 
Observations collected until 2003 were consistent with a sce- 
nario in which the twist angle was steadily increasing before 
the glitch epochs, culminating with glitches and a period of in- 
creased timing noise, and then decreasing, leading to a smaller 
flux and a softer spectrum. Both the Chandra and the Swift ob- 
servations caught the source in a (relatively) hard, luminous 
state, revealing a reversed trend. However, the hardening-flux 
correlation is maintained, lending further support to this sce- 
nario. 

What is particularly interesting, and measured here for the 
first time, is that since the last XMM-Newton observation the 
fraction of the total flux in the power-law component slightly 
decreased although the source spectrum became harder. This 
is somehow counter-intuitive. If taken at face value, it may 
be explained by the fact that, in the twisted magnetosphere 
model, both the spatial distributions of the magnetospheric cur- 
rents (which act as "scattering medium") and the surface emis- 
sion induced by the returning currents (which acts as source 
of seed photons for the resonant scattering) are substantially 
anisotropic. Seed thermal photons and scatterers are confined 
in two different limited ranges of magnetic colatitudes, and 
both distributions move away from the poles for larger twist 
angle, although at a different rate. For instance, by using the ex- 
pressions provided by Thompson, Lyutikov & Kulkarni (2002) 
for the differential luminosity induced by the returning cur- 
rents, we can estimate that the center of the heated surface re- 
gion moves from ^ 37° to ~ 63° in colatitude when A0ns in- 
creases from ~ 0.1 to 2 radians. Correspondingly, the peak of 
the efficiency of the scattering only shifts from ^ 66° to ~ 72° 
in colatitude. The size of the region interested by the scatter- 
ing decreases to ^ 37%, while the thermally emitting region 
becomes ^ 5% larger. Although the model is quite approxi- 
mated, and the above numbers should be treated with care, this 
strong anisotropy suggests that the observed drop in the non- 
thermal flux may be due to the fact that a lower fraction of soft 
photons are intercepted by the cloud of scattering particles sur- 
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Fig. 2. Left Panel: Chandra unfolded spectrum fitted with an absorbed blackbody plus a power-law. See Table 2 for further details 
on the spectral parameters. Right Panel: Chandra pulse profile in the 0.3-8.0keV energy band. 
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Fig. 3. Long term spectral history of IRXS J170849. 0^00910, correlated intensity-hardness variation (adapted from Rea et al. 
2005). All reported fluxes are unabsorbed and in the 0.5-lOkeV energy range. For clarity, the observations dates are: ROSAT- 



1994; ASCA - 1996; first BeppoSAX- 1997; second B eppoS AX - 
CC- 2004; Swift -2005. 



2001 ; Chandra-HETG - 2002; XMM-Newton - 2003; Chandra- 



rounding the star. Clearly, since the scattering depth increases 
with the power-law will be in any case flatter. 

Our preliminary quantitative estimates show, however, that 
the increase in size of the thermally emitting region is not suf- 
ficient to account for the observed variation in the blackbody 
radius, at least on the basis of the original model by Thompson, 
Lyutikov & Kulkarni (2002). We only note in this respect that 
viewing geometry effects may be important, since the expected 
change in the position of the heated surface region may result 
in a larger portion of the emitting area coming into view. 

Finally, we might speculate that the long-term variations 
shown in Fig. [3] may have a cyclic behavior with a recurrence 
time of « 5-10 yr. A possible explanation within the magnetar 
scenario might be the periodic twisting/untwisting of the star 
magnetosphere, where the characteristic dissipation time of a 
static twist is in fact « 1-10 yr according to the more recent es- 
timates (Beloborodov & Thompson 2006). A detailed study of 
this intensity-hardness correlation, through further X-ray mon- 



itoring of this source, is needed in order to better constrain the 
model, and to infer information on the physical conditions in 
the star magnetosphere. Note that with a detailed modeling of 
this correlation we would be able in the near future to predict 
the occurrence of glitches and possibly also of bursts. 
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